Cosmological Parameters from Pre-Planck CMB Measurements 
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Recent data from the WMAP, ACT and SPT experiments provide precise measurements of the 
cosmic microwave background temperature power spectrum over a wide range of angular scales. 
The combination of these observations is well fit by the standard, spatially flat ACDM cosmological 
model, constraining six free parameters to within a few percent. The scalar spectral index, n s = 
0.9690 ± 0.0089, is less than unity at the 3.5o level, consistent with simple models of inflation. The 
damping tail of the power spectrum at high resolution, combined with the amplitude of gravitational 
lensing measured by ACT and SPT, constrains the effective number of relativistic species to be 
iVeff = 3.28 ± 0.40, in agreement with the standard model's three species of light neutrinos. 



Introduction — It has long been appreciated that the 
Cosmic Microwave Background (CMB) power spectrum 
contains enough information to precisely determine the 
standard model of cosmology [IH1]. This promise has 
been realized through a series of increasingly sensitive 
experiments, most recently with the WMAP satellite's 
9-year full-sky observations [6] and the arcminute- 
resolution maps from the Atacama Cosmology Telescope 
(ACT) [7HS] and the South Pole Telescope (SPT) [101 [IT]. 
The combination of these measurements probes the tem- 
perature power spectrum on angular scales ranging from 
90 degrees to 4 arcminutes, scales at which the primary 
cosmological temperature fluctuations dominate. The 
primordial fluctuations are well approximated as a Gaus- 
sian random field E2], but ACT and SPT have also 
detected the non-Gaussian features due to gravitational 
lensing of the microwave radiation by the intervening 
large-scale structures |T3l [14] . In this letter we present a 



joint analysis of the ACT, SPT, and final WMAP 9-year 
power spectra to obtain an estimate of the cosmological 
parameters from microwave background data alone. 

Data and Analysis Method — In Figure [T] we show the 
compilation of CMB temperature power spectra used in 
this analysis. At large angular scales we use the temper- 
ature and polarization data, and associated likelihood, 
from the 9-year WMAP analysis [hereafter WMAP9, 
[5]. This measures the Sachs- Wolfe plateau and the first 
three acoustic peaks, 2 < I < 1000. At smaller scales, 
500 < I < 3500, we use data from ACT and SPT. 

Here we follow the method introduced in [8] to estimate 
the primary CMB bandpowers from both sets of spectra, 
marginalizing over the possible additional power from 
Galactic and extragalactic emission, and the Sunyaev- 
Zel'dovich effects. We use a Gibbs sampling method to 
simultaneously estimate CMB bandpowers and a set of 
ten secondary parameters. For ACT we extract primary 
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FIG. 1. WMAP9 temperature data and ACT and SPT CMB lensed bandpowers marginalized over secondary emissions. The 
ACT bandpowers are estimated separately for ACT-S and ACT-E and coadded here with an inverse variance weighting. The 
SPT bins are highly correlated, (50 — 65%) at small scales, I > 2000, due to foreground uncertainty. The correlation is about 
5% between neighbouring ACT bins. The solid line shows the lensed CMB best fit obtained combining the three datasets. The 
ACT and SPT bandpowers are available on LAMBDA | |http:// lambda.gsfc.nasa.gov/|). 



CMB bandpowers from the 148 and 218 GHz auto and 
cross power spectra from two regions (ACT-E and ACT- 
S, [I]) °f the sky [T5], taking the multi- frequency band- 
powers in the range 500< I < 10000. We include SPT 
150 GHz data [10] from 650< I <3000, and marginalize 
over a common model for secondary components |16j . We 
impose a Gaussian prior of 12. 3 ±3. 5 /iK 2 at i = 3000 on 
the SPT radio source Poisson power, having subtracted 
7 /xK 2 of cosmic infrared background Poisson power, 
treated separately in our likelihood, from the total ex- 
pected Poisson level [10], [17] • The resulting ACT and 
SPT lensed bandpowers are shown in Figure [T] and the 
secondary parameters are consistent with those reported 
in [5] [9] . The errors shown are the diagonal elements of 
the covariance matrix, with the SPT calibration error re- 
moved for consistency with ACT. The full covariance ma- 
trix includes correlations due to foreground uncertainty, 
beam error, and the overall calibration for SPT. 

We then construct an ACT+SPT likelihood from these 
CMB bandpowers, which can also be used for each ex- 
periment on its own. This is a Gaussian distribution us- 
ing 42 data points from ACT (21 each from ACT-E and 
ACT-S) and 47 from SPT, with an associated covariance 
matrix. For ACT we only use I <3500 bandpowers in the 
likelihood, where their distributions are Gaussian. When 
combining ACT with SPT, we use only ACT-E data to 
eliminate the covariance between ACT-S and SPT, which 
observe overlapping sky regions. We combine this like- 
lihood with WMAP9, using the CosmoMC code [T5] to 
estimate cosmological parameters. 



TABLE I. Standard ACDM parameters from the combination 
of WMAP9, ACT and SPT. 



Parameter 


WMAP9 


WMAP9 


WMAP9 




+ACT 


±SPT 


±ACT±SPT a 


100^6^ 


2.260 ±0.041 


2.231 ±0.034 


2.252 ±0.033 


l00Q c h 2 


11.46 ± 0.43 


11.16 ±0.36 


11.22 ±0.36 


1006U 


1.0396 ± 0.0019 


1.0422 ±0.0010 


1.0424 ±0.0010 


T 


0.090 ±0.014 


0.082 ±0.013 


0.085 ±0.013 


n s 


0.973 ±0.011 


0.9650 ± 0.0093 0.9690 ± 0.0089 




2.22 ±0.10 


2.15 ±0.10 


2.17 ±0.10 


a A " 


0.716 ±0.024 


0.737 ±0.019 


0.735 ±0.019 




0.830 ±0.021 


0.808 ± 0.018 


0.814 ±0.018 


to 


13.752 ±0.096 


13.686 ±0.065 


13.665 ±0.063 


H 


69.7 ±2.0 


71.5 ± 1.7 


71.4 ± 1.6 


lOOr-s/Dvo.57 


7.50 ±0.17 


7.65 ±0.14 


7.66 ±0.14 


100r./£Vo.35 


11.29 ± 0.31 


11.56 ±0.26 


11.57 ±0.26 


best fit X* 


7596.0 


7617.1 


7640.7 



a The combination ACT+SPT uses ACT-E data only. 
We report errors at 68% confidence levels. 

b Derived parameters: Dark energy density, the amplitude of 
matter fluctuations on 8 ?t _1 Mpc scales, the age of the 
Universe in Gyr, the Hubble constant in units of km/s/Mpc, 
and the galaxy correlation scales at redshifts 0.57 and 0.35. 

We consider the basic spatially-flat ACDM cosmologi- 
cal model defined by six parameters: the baryon and cold 
dark matter densities, f2f,/i 2 and f2 c /i 2 ; the angular scale 
of the acoustic horizon at decoupling, 9a] the reionization 
optical depth, r; the amplitude and the scalar spectral 
index of primordial adiabatic density perturbations, A 2 ^ 
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FIG. 2. Marginalized one-dimensional distributions for the 
six basic ACDM parameters, for combinations of WMAP9 
(W9), ACT (A) and SPT (S) data. 



and n s (at a pivot scale fco = 0.05 Mpc~ x ). We also ex- 
tend the standard model to include a seventh parameter 
iV e ff , the effective number of relativistic species at decou- 
pling. The high-£ damping tail measured by ACT and 
SPT is particularly sensitive to this parameter. 

Results and discussion — The simple ACDM model fits 
all the data well, with the estimated parameters shown in 
Tableland Figure § We find x 2 /dof= 37.9/42 (Prob- 
ability to Exceed, PTE = 0.65) for ACT and 53.2/47 
(PTE = 0.25) for SPT when combined individually with 
WMAP9, assuming the degrees of freedom equal the 
number of additional data points. The best-fitting pa- 
rameters for ACT are all within about la of the corre- 
sponding best-fitting parameters for SPT. For the com- 
bined analysis, the ACT+SPT best fit x 2 / d <ff is 78.9/68 
(PTE = 0.17). Compared to the joint best-fit model, the 
SPT-only best-fit has A X 2 = 2.5 worse and the ACT- 
only best fit (for ACT-S+ACT-E) has A X 2 = 2.2, indi- 
cating that the common model fits both datasets. Figure 
[3] shows the residual power for the high-€ datasets after 
subtracting the joint best-fitting model. We do not ob- 
serve any particular features in ACT; the SPT power is 
more suppressed at multipoles I > 1500, but the points 
include an uncertain correlated extragalactic foreground 
contribution, whose dominant term is a Poisson shape. 

The addition of ACT and SPT helps WMAP constrain 
the basic six parameters due to a more precise determi- 
nation of the higher order acoustic peak positions and 
amplitudes. The measurement of 9a improves by a fac- 
tor 2.2 and the error on the baryon density is 1.6 smaller 
compared to WMAP9 alone. However, as noted in [TT] . 
the increased acoustic horizon scale leads to a predicted 
distance, Dy, to objects at redshift z — 0.57, in units of 
the sound horizon at recombination, r s , of I00r s /Dy — 
7.66 ± 0.14, more than 2a larger than measured by the 
BOSS experiment ([15], 100r s /D v = 7.3 ± 0.1). The 
prediction at z = 0.35, 100r s /D v = 11.57 ± 0.26, is 
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FIG. 3. Residual power after subtracting the same best-fitting 
lensed CMB model. The reduced X 2 / dof for ACT is 40.1/42 
(PTE=0.55) and for SPT 55.7/47 (PTE=0.18). We show 
ACT-E and ACT-S coadded residuals. The grey band in the 
bottom panel shows the 2a uncertainty in the Poisson source 
component. Overall calibration errors are not included. 



consistent at la with the SDSS DR-7 observations f |20j. 
100r s /Dv — 11.3±0.2). We find a preference for a scale- 
dependent primordial power spectrum at 3.5c from the 
CMB, with n s = 0.9690 ± 0.0089 at 68% confidence. 

The CMB power spectrum is sensitive to the compo- 
sition of the Universe. The radiation energy density is 
the energy density in photons plus the sum of the energy 
density in relativistic species that do not couple electro- 
magnetically, including standard model neutrinos. We 
parametrize the energy density in other relativistic par- 
ticles through N c g. In the standard cosmological model, 
N e g = 3.046 [2~lTf2"3"] describes the three known neutrino 
species. If there is an extra neutrino species that decou- 
ples at the same temperature as the standard neutrinos 
then 7V e ff — 4. If, instead, there is another light weakly 
interacting stable particle that decouples earlier, it will 
increase N e g by the cube of the ratio of the decoupling 
temperatures. The extra energy density in relativistic 
species has three noticeable effects on the CMB power 
spectrum [24H27] : (1) it increases the expansion rate of 
the Universe, which impacts both the acoustic and damp- 
ing scale, an effect that is mostly degenerate with increas- 
ing the matter density, 0, m h 2 ; (2) it modulates the helium 
abundance from big bang nucleosynthesis, which in turn 
modifies the damping tail through free electrons avail- 
able at recombination; and (3) the relativistic particles 
will free stream out of density fluctuations and suppress 
the amplitude of the power spectrum on small angular 
scales (an effect partially degenerate with increasing n s , 
which enhances the amplitude of the power spectrum on 
those scales). 

In our analysis, we vary N c g as a free parameter and 
we assume that the same relativistic species are present 
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FIG. 4. Left: Marginalized distribution of N B s for different data combinations, showing consistency with three neutrino species. 
Middle and Right panels: Marginalized 68% and 95% contours in the N c a- n s and N c a- Ho planes; Neff is correlated with both 
parameters. The standard model expectation of N e g = 3.046 is indicated with dashed lines. 



at nucleosynthesis. We keep other quantities that de- 
scribe the damping tail set to standard values: the total 
available electron abundance (determined by the primor- 
dial helium abundance Y p ) is consistent with standard 
big bang nucleosynthesis; the higher precision determi- 
nation of recombination is used [2"5 1, as implemented in 
[29| . Combining the data we are considering in this work 
we find at 68% confidence level: 

iVeff = 2.90 ± 0.53 (WMAP9 + ACT) 
Nes = 3.75 ± 0.47 (WMAP9 + SPT) 
N eS = 3.37 ±0.42 (WMAP9 + ACT + SPT). 

In Figure [4] we show the distribution for N e g from 
WMAP9 combined with ACT and SPT separately, and 
together. There is a 1.2er difference between the ACT 
and SPT estimates; as noted in [IT], the SPT data prefer 
a higher value, indicating more suppression of the small- 
scale spectrum. The probability that this variation is 
given by statistical scatter is around 50%, [3D]. Based 
on the difference in the damping tail measurements, [32] 
decided to not combine the ACT and SPT data. In this 
paper we take a different approach and view the con- 
sistency sufficient for combination. As noted in [33], a 
Bayesian model comparison shows no evidence in favor 
of adding an additional N e g parameter beyond those of 
the standard cosmology. 

The correlations among N e g, n s , and Hq are also 
shown in Figure [4] the suppression of small scale power 
due to larger values of N e g can be partially compensated 
by increasing n s and tt m h 2 . This leads to a larger derived 
value of Hq if the CMB peak positions are held fixed. 

Since a higher value for N e g requires a higher mat- 
ter density today, it gives a higher amplitude of grav- 
itational potential fluctuations and an increased gravi- 
tational lensing signal. Measurements of the four-point 
function of the CMB temperature maps provide a mea- 
surement of the lensing deflection signal. The ACT [7] 
and SPT [14] data constrain the amplitude of the lensing 



potential power spectrum at t = 400 to be C£q — (3.69± 
0.80) x 10~ 8 for ACT and C^ = (2.92 ±0.54) x 10~ 8 for 
SPT, yielding a combined C| K = (3.17 ± 0.45) x 10~ 8 . 
Adding this constraint gives 



N, 



eff 



3.28 ± 0.40 (WMAP9 + ACT + SPT + Lensing), 



consistent with three neutrino species. 

Conclusions — Current microwave background power 
spectrum measurements are consistent with the standard 
ACDM cosmological model, and independent data sets 
are consistent with each other, with a mild tension be- 
tween the ACT and SPT damping tails. Upcoming maps 
from the Planck satellite will provide independent mea- 
surements of the same sky regions with excellent fore- 
ground characterization. 
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